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^ ' Abstract. We discuss the three-body decay mechanisms of many-body reso- 

^ , nances. R-matrix sequential description is compared with full Faddeev com- 

putation. The role of the angular momentum and boson symmetries is also 
studied. As an illustration we show the computed a-particle energy distribu- 
^ . tion after the decay of "'^^0(1"^) resonance at 12.7 MeV. 
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The three-body decay of many-body resonances can be accurately measured in 
complete kinematics. Information about the decaying state and the decay mech- 
anism is usually extracted from the measurement of the three fragments after 
the decay. Although this is a common practice, the situation is ambiguous. The 
^ ' experimental analyses of these processes are based on the i2-matrix formalism 

^ . which inherently assumes two successive two-body decays. The input are the 

properties of the intermediate two-body states and the population of the nuclear 
many-body initial state approximated as a three-body system. 

Occasionally, in principle contrary descriptions are able to explain the ob- 
served distributions making the understanding of the underlying physics diffi- 
cult. An example demonstrating the difficulties is the 3a decay of l"*" state in 
^^C which was successfully described by two opposite mechanisms: a sequential 
decay via the state in ^Be [1], and a direct decay into the three-body con- 
tinuum [2]. This requires an explanation. What information is contained in a 
full-kinematics measurement of three-body decay? Apparently unique informa- 
tion can only be extracted under favorable conditions. The crux of the matter is 
that the decay mechanism is related to a "decay path" , an intermediate structure, 
which in contrast to the final state signal in the detector is not an observable. 
We shall in this contribution compare the results from three-body calculations 
and experimental i?-matrix analyses. 
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Three-body decays: structure, decay mechanism and fragment properties 



2 Energy distributions 



The large-distance observable structure of the many-body initial state is a three- 
body continuum state, therefore we compute the resonance structure in a three- 
body cluster model [3]. We use the complex scaled hyperspherical adiabatic ex- 
pansion method to solve the Faddeev equations which describe the 3-body sys- 
tem. The appropriate coordinates are the so-called hyperspherical coordinates 
and consist of the hyperradius = 4^^^^(rj — RY , and five hyperangles. In 
the adiabatic hyperspherical expansion method the angular part of the Faddeev 
equations is solved first and the angular eigenfunctions 'PnJM are then used as a 
basis to expand the total wave function ^■^^^ . 

We include short-range [4] and Coulomb potentials. The many-body effects 
that are present at short distances are assumed to be unimportant except for 
the resonance energy. This is taken into account by using a structureless 3-body 
interaction that fits the position of the resonance. 

The resonance wave-function contains information about the decay mecha- 
nism, and the large-distance properties reflect directly the measurable fragment 
momentum distributions. The single particle probability distributions are ob- 
tained after integration of the absolute square of the wave function over the four 
hyperangles describing the directions of the momenta. 

The many-body initial state res- 
onance evolves into three clusters at 
large distances. The total angular mo- 
mentum and parity J'^ is conserved 
in the process. This symmetry com- 
bined with Bose-statistics imposes con- 
straints on the resulting momentum 
distributions. An early example of these 
effects applied to three pion decays can 
be found in ref. [5]. Fig. [1] shows the re- 
gions of the Dalitz plot where the den- 
sity must vanish for the decay of a l"*" 
state into three a-particles. This gives 
rise to the minima in the single a-particle energy distribution. The Dalitz plot 
computed within the Faddeev framework is also shown and is in agreement with 
these symmetry constraints. 

Fig. [2] shows the single-a energy distributions, i.e. the probability for emer- 
gence of one a particle with a given energy divided by its maximum allowed, for 
the 1^ state of ^^C computed with i?-matrix analysis [1,6]. This corresponds to 
the projection of the Dalitz plot in fig.[T]on the y-axis. The decay is assumed to be 
sequential via *Be(0"'") since angular momentum forbids the decay via ®Be(2"'"). 
We have varied the two-body energy and width. When both the two-body energy 
and width are small (left) a narrow peak corresponding to the emission of the first 
a arises. The other two a's are related to the broad peaks. By increasing the two- 
body energy and width the three-peak distribution becomes rather pronounced 
and insensitive to the two-body parameters when either E2r/E^r is larger than 
about 0.5 or the two-body width is large. The same figure contains the curves 




Figure 1. (Color online) Regions of the 3q 
Dalitz plot where the density must vanish (in 
black) for a 1"*" state (left). Dalitz plot for 
the 1"*" state of ^^C. x-axis corresponds to 
(_BQi/2£a2)/\/(3) and y-axis to E^i in MeV. 




Figure 2. (Color online) Single q energy distributions from R-matrix analysis for the decay 
into three a-particles of the ^■^C(l^) resonance at 12.7 MeV of excitation energy. The energies 
and widths of the intermediate *Be(2"'") state are varied as specified in the panels. The green 
dashed curve corresponds to the case where the symmetrization of the wave function is omitted. 
The middle panel corresponds to the measured resonance energy. 



Table 1. The probability Peeq for populating the component related to the decay via *Be(2"'") 
at large distances in the computation of the 1"*" resonance of ^^C for a complex rotation angle 
9 = 0.25. The three-body energy (£3, ) is varied by adjusting the strength of the three-body 
potential. The *Be(2"'") two-body energy is maintained E2r ~ 2.7 MeV. The energies are referred 
to the 3a or 2a separation threshold. 
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corresponding to the case where the boson symmetry is omitted. For a low and 
narrow two-body state the effect of this symmetry seems to be unimportant, but 
an increase on the width leads to a two-peak (not three-peak) distribution. 

Fig.[3]shows the results from the full three-body computation and the compu- 
tation from the lowest continuum three-body wave function (K=8) from ref. [7] 
(democratic decay). This is the simplest assumption with the correct symme- 
tries. We have varied the three-body energy and consequently the three-body 
width. Two rotation angles have been considered: one of them is large enough 
to accumulate the contribution of sequential decay through the ®Be resonance in 
a single adiabatic potential, while the other is not. The calculations include the 
boson symmetry of the a-particles. The results from the large rotation angle do 
not include the contribution from the decay via ^Be(2"'') and are very close to 
the democratic decay. In the result from the full Faddeev computation the three 
peaks are closer to each other and this approaches better the experiment. The 
fractions of population at large distance are given in table [2] for the different val- 
ues of £^3r- shown in fig. [3l We can observe that the sequential decay probability 
increases as we increase the three-body energy. 




Figure 3. (Color online) Single a energy distributions from Faddeev computation for direct 
decay of the ^^0(1+) resonance. The three-body energy is varied by changing the strength 
of the three-body potential. The relative energies and widths of the intermediate *Be(2"'") 
resonance are specified in the panels. The solid (black) and dashed (green) curves correspond 
to complex rotation angles of ^ = 0.25 and 0.1 respectively. The dotted (blue) curve corresponds 
to democratic decay [7]. For Q — 0.25 only direct decay is shown. 

3 Conclusions 

We have computed the observable momentum distributions from decay of three- 
body resonances by use of i?-matrix simulations and from full Faddeev calcula- 
tions. We have considered the example of the l"*" resonance of ^^C. The angular 
momentum and boson symmetries constrain the resulting momentum distribu- 
tions. The same measured momentum distributions can be described in differ- 
ent complete basis sets, e.g. either direct products of two-body states and their 
center-of-mass motion relative to the third particle (i?-matrix) or three-body 
continuum wave functions (Faddeev). The fact that different descriptions seem 
to work indicates that the same wave function could be described in different 
ways. Extracting information of both structure and decay mechanism can then 
be misleading and requires model interpretations. Full Faddeev computations 
successfully reproduce the measured distributions. 
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